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Queuing models provide insight into the temporal inhomogeneity of human dynamics, charac- 
terized by the broad distribution of waiting times of individuals performing tasks. We study the 
queuing model of an agent trying to execute a task of interest, the priority of which may vary with 
time due to the agent's "state of mind." However, its execution is disrupted by other tasks of random 
priorities. By considering the priority of the task of interest either decreasing or increasing alge- 
braically in time, we analytically obtain and numerically confirm the bimodal and unimodal waiting 
time distributions with power-law decaying tails, respectively. These results are also compared to 
the updating time distribution of papers in the arXiv.org and the processing time distribution of 
papers in Physical Review journals. Our analysis helps to understand human task execution in a 
more realistic scenario. 

PACS numbers: 89.20.-a, 89.75.Da, 89.65.-s 



I. INTRODUCTION 

Understanding human behavior is a fascinating and 
challenging subject that has been explored extensively 
in the fields of sociology, economics, and psychology for 
centuries [IH3]. Recently, with the need to gain deeper 
insight into complex social systems, the human behavior 
has also been studied from the perspective of other scien- 
tific disciplines including physical and computational sci- 
ences [Ij. This is because the concepts and methods de- 
veloped for investigating physical processes have turned 
out to be applicable to understand the human and social 
phenomena. Moreover, the huge amounts of digital data 
have enabled us to explore the human behavior at the 
very high resolution. The empirical studies based on such 
dataset have shown that the timing of human actions is 
not simply random but correlated or bursty f&H^. Such 
correlations are often characterized by the heavy tailed 
or power-law distribution of waiting or response time r as 
P(r) ~ T~", where awl for e-mail communication and 
1.5 for letter correspondence, respectively OIH]. Here r 
is defined as the time taken by the user to respond to a 
received message. 

The origin of bursts in human dynamics has been ex- 
plored for years. The long inactive periods at nighttime 
as well as on weekends are known to result in the bursts 
of e-mail communication [51 llOj , where the circadian and 
weekly cyclic patterns of humans were modeled by the 
non-homogeneous Poisson process. In addition to these 
human cyclic patterns, the human dynamics can be also 
affected by other human factors, such as the task execu- 
tion behavior. To figure out the role of task execution 
behavior we have adopted the framework of queuing the- 
ory [SI [TTHTB] , where the to-do list of an agent is modeled 
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as a finite length queue. The agent assigns priorities to 
tasks in the queue and then executes them according to 
some protocol. As in the Barabasi's model [5j, if the task 
with highest priority is executed first, then the one with 
low priority must wait for a relatively long time in the 
queue, leading to the heavy tailed distribution of waiting 
times. Barabasi's priority queuing model is based on a 
few key assumptions: (i) It is a single agent model, where 
the agent has no explicit interaction with other agents, 
(ii) all the tasks in the queue are of the same type, such 
as replying to the received message, and (iii) the priority 
assigned to each task is kept fixed over time. Several vari- 
ants of this model have been studied by means of relaxing 
at least one of these assumptions [T^17) . Here we per- 
form the similar study by relaxing the assumptions (ii) 
and (iii) and consider that only one type of the task has 
the time- varying priority. In the similar context, Blan- 
chard and Hongler studied the queuing systems on the 
basis of population dynamics, where the individuals in 
a city (tasks in a queue) are assigned with the priorities 
increasing with time due to aging or the deadline effect, 
and finally dying (executed) according to "the highest 
priority first" protocol [14]. However, instead of the pop- 
ulation dynamics based models, the microscopic models 
with simpler setups might help us to better understand 
the effect of time-varying priorities on the waiting time 
distributions. 

In this paper, we adopt the microscopic approach and 
study a single agent queuing model with two tasks: one 
is the task of interest with time- varying priority and the 
other is the unexpected and random task with random 
priority. As an example let us consider a situation in 
which a researcher wants to complete the research project 
of current interest, but he or she gets continuously dis- 
tracted by meetings, coffee breaks and so on. The priority 
of the task of interest recognized by the researcher can 
vary due to his or her state of mind, such as mood and 
happiness. The priority may increase in some cases, such 
as when the deadline assigned to the task approaches [18] , 
or decrease in other cases, for example, when the re- 
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FIG. 1: Schematic diagram of time- varying priority queuing 
model, (a) The task of interest (type-A) has the time- varying 
priority x{t), which is compared with the random task (type- 
B) with random number rt at each time step and executed if 
x{t) > rt- (b) The case with power-law decreasing priority, 
(c) The case with power-law increasing priority: bounded and 
unbounded cases. See the text for details. 



searcher gradually loses interest in that task. For conve- 
nience, we call the task of interest and the random task as 
type-A and type-B tasks, respectively. This model with 
one type-A and one type-B tasks is minimal and can be 
generalized to more complex scenarios, such as the case 
of an agent with many tasks of various types or the case 
with interacting agents. 

The paper is organized as follows. In Section |llj we 
study the models with different cases of time- varying pri- 
orities. For each case we obtain the analytic solutions for 
the waiting time distribution that are confirmed by nu- 
merical simulations. In Section |III[ we show how the 
extended versions of our models to many agents can be 
used to understand the paper updating mechanism in 
arXiv.org and the paper reviewing process in Physical 
Review journals. Finally we summarize the results and 
make conclusions in Section Hvl 



II. MODEL 

We consider an agent with a queue of size 2 as depicted 
in Fig. TJa). The first site of the queue is occupied by 
the task of interest with time-varying priority (type-A 
task). When the type-A task is introduced to the queue 
at the time step t = 0, its priority is given by the ran- 
dom number Xq drawn from the uniform distribution of 
(0, 1) and then it varies with time. The second site of 
the queue is for the random task (type-B task), which is 
replaced by a new random task at every time step. The 
priority of the type-B task introduced at time step t is 
given by the random number rt drawn from the uniform 
distribution of (0,1). At each time step i, the priority 
x{t) is compared with the random number rt and exe- 
cuted when x{t) > rt- Here the execution time defines 
the waiting time r for which the type-A task waits for 



the execution since it is introduced to the queue. This 
is the deterministic case where only the highest priority 
task is executed. 

We also consider a stochastic version of the model, 
where the lower priority task is allowed to be executed 
with small probability. At each time step, the highest pri- 
ority task is executed with the probability p, while with 
probability 1 — p one of two tasks is selected at random 
for the execution, which we call random selection. This 
implies that the priority is not the only determinant for 
the execution. This can also be interpreted as the "trem- 
bling hand" effect in the game theory (TH]. In general, 
the waiting time distribution P(t) is obtained by 
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where e = is assumed to be very small in this paper. 



A. Fixed priority 

We first consider the case of fixed priority. Here the 
priority of the type-A task is fixed, i.e. x(t) = xq Vi. The 
waiting time distribution is given by 



P(t) = p- 



T+1 



In the limit of r ^ 1, it can be approximated by 



P(r) « (r-2 + 2er-i) 
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where we define the cutoff of distribution as Tc{p) — 
[ln(l/p)]"i « (2e)-^ In the scaling regime, i.e. 1 <C 
T ^ Tc ~ (2e)~^, the first term r^^ dominates the distri- 
bution function, so that P{t) ss r^" with a = 2. In the 
deterministic case with p — I, P{t) ~ {t+i){t+2) ~ "^^^ 
in the limit of t ^ 1. Note that the small probability of 
random selection does not change the scaling exponent 
but leads to the finite cutoff of the distribution. Figure [2] 
shows that the analytic solution is confirmed by the nu- 
merical simulations for both deterministic (p = 1) and 
stochastic {p = 0.999) cases. 



B. Power-law decreasing priority 

Now we consider the case when the priority of the type- 
A task changes over time. This change might be caused 
by the internal factors related to the agent as well as to 
the external factors. We consider the simple case with 
the algebraically decreasing priority as 
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FIG. 2: The case with fixed priority as x{t) = xq for deter- 
ministic (p — 1) and stochastic {p = 0.999) versions. For the 
numerical results we used 10® samples. The solid line for the 
stochastic version represents the analytic solution of Eq. 



The decreasing speed of the priority is controlled by the 
exponent 7. If 7 = 0, the model reduces to the fixed 
priority case, where the waiting time distribution decays 
as the power-law with exponent a = 2. For the larger 
value of 7, the priority decays faster so that the type-A 
task becomes less likely to be executed. Hence one may 
expect that this will lead to a heavier tail of waiting time 
distribution and consequently to a smaller value of a. 

At first we analyze the case with 7 <C 1. Since the 
priority x{t) decreases very slowly, we can use the ap- 
proximation of x{t) « x{t) — j^^r^Tiyf- Putting this into 

Eq. (g, we get 
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FIG. 3: The case with power-law decreasing priority as 
^(^) ~ (t+ip ^'^^ deterministic version with p = 1 (a) and 
for stochastic version with p = 0.999 (b). For the numerical 
results we used up to 10* samples. In the inset of (a) we plot 
the numerical results of 0(7) by the red filled circles with the 
analytic expectation by the black lines. In (b) the solid lines 
represent the analytic solutions of Eq. Q for 7 — 0.2 and 
Eq. ||9| for 7 = 2, respectively. 
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For large values of r and with the fact that is close to 
1, the above equation reduces to 



(6) 



In the scaling regime of 1 ^ r ^ w (2e)~^, the first 
term t"'^^"''' dominates the distribution function, so that 
P(r) w T^" with a = 2 — 7. Thus the exponent a is 
smaller for larger value of 7, as expected. 

For the deterministic version of the model, i.e. p = 1, 
we perform the numerical simulations with various values 
of 7 as shown in Fig. |3]^a). We find that for any value of 
7 the waiting time distribution follows the power-law be- 
havior with the power-law exponent a as a function of 7, 
see the inset of Fig.jsjja). Our analytic result of a = 2 — 7 
is numerically confirmed up to 7 ~ 0.6. The numerical 
results, however, suggest that the value of a has the min- 
imum at 1.22 for 7 = 1, while it increases when 7 > 1. In 



the stochastic case with p < 1, the second term r^'^^''') 
in Eq. ([6| begins to affect the dominant power-law behav- 
ior and the exponential cutoff also becomes finite. This 
approximate solution is again numerically confirmed for 
various values of 7, as depicted in Fig. [sjjb). 

Next, the other limiting case when 7 ^ 1 is considered. 
Since the value of x{t) decreases very fast, we use the 
approximation of x(t) <^ 1 for all t. Thus the product 
term in Eq. ([I]) can be expanded in the limit of e ^ 1 as 
follows 
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Then the waiting time distribution is obtained by 
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Using the assumption of 7 ^ 1 and the condition that 
er <C 1, the above equation can be approximated as 



P(r) 
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(10) 



For the scaUng regime, the waiting time distribution is 
dominated by the first term, leading to P{t) « r^" with 
(3 = 7. This is numerically confirmed up to 7 > 1.4 in 
the deterministic case as shown in the inset of Fig. |3][a). 
For the stochastic case with p = 0.999, we find that the 
numerical simulations confirm the analytic solutions but 
with some deviation due to the approximation we have 
used, see Fig. [sjjb). 



The result of Eq. ( 10 1 can be understood more intu- 



itively as following: For large r, the probability that the 
type- A task is not executed becomes of the order of 1, 
i.e. nr=o^[l~^(0 + ^] ~ ^(1)- Therefore, the only factor 
relevant to determine P{t) is the probability of the task 
being executed at time step r. Thus, 



P(r) ^ [x(t) + e]e-^/^^ « [t'T + e]e~^''% 



(11) 



which implies = 7. For the stochastic case we find 
that the second term in the result is only a function of e 
not coupled to r. That is, the waiting time distribution 
does not follow a typical power-law with an exponential 
cutoff as r^"e^'^/'^'=, but shows a bimodal combination 
of the power-law and exponential distributions. Similar 
bimodal distribution has been empirically observed for 
the inter-event time distributions of Short Messages in 
mobile phone communication [5] . 



C. Power-law increasing priority 

In this Subsection, we consider the case of the priority 
of the task increasing as a power-law with time, starting 
from a;(0) = 0. The priority can increase up to some 
bounded value, which we call the bounded case. Here the 
priority of type-A task is given by 
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where the scale factor d plays the role of a deadline given 
to the task, xq can be interpreted as an intrinsic prior- 
ity of the task and its value is drawn from the uniform 
distribution of (0, 1). 

On the other hand, the priority of the type-A task can 
increase up to the maximum possible value of 1, which 
we call the unbounded case. Once the priority reaches 1, 
the task is inevitably executed. In this case, the priority 
x{t) is given by 
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1 \it>tr. 



(13) 



where tr 



d is defined by the condition x{tc) — 1. 
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FIG. 4: The case with power-law increasing priority as 
x(t) = a;o(|)'' for deterministic version with p = 1 (a) and 
for stochastic version with p — 0.999 (b). The value of x{t) 
can increase to xq in the bounded case, where x{t) = xo for 
t > d. On the other hand the value of x{t) can increase to 
1 in the unbounded case. For the numerical results we used 
up to 10^ samples. In (b) the solid lines for 13 — 2 represent 
the analytic solutions of Eq. ( [T4| | for r < 30 and Eq. ( [T6| for 
T > 150, respectively. 



For both the bounded and unbounded cases, we can 
adopt the approximation used in Eq. ([7| in the limit of 
r <C d, resulting in 
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In the limit of er ^ 1, the above equation can be ap- 
proximated as 



(15) 



For the range of e^^^ <^ ^ <^ 1, the waiting time distri- 
bution is dominated by the first term, indicating a power- 
law increasing behavior as P{t) ~ r" with a' = /?. This 
solution is compared to the numerical results for both the 
deterministic and stochastic versions of the model, where 
we set d = 100, see Fig. |4| 
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On the other hand, in the hmit of r ^ d, the waiting 
time distributions decay as a power-law, but with dif- 
ferent tail exponents for bounded and unbounded cases. 
In the bounded case, when t > d, the priority is fixed 
as x(t) = xq, corresponding to the fixed priority model. 
Therefore, we get a — 2, which is independent of the 
values of /? and d. To obtain the distribution for the un- 
bounded case, we use the approximations that x{t) <^ 1 
if t < d and that x{t) ~ x(t) if t > d. In addition, the 
type-A task always gets executed as the priority reaches 
1, implying that x{t) + e < 1. This means that the value 
of xo has an upper bound of Xc given by Xc{2)^ + e — 1- 
Thus the upper limit of the integral range in Eq. ([T]) is 
also limited to Xc, leading to 
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The above 



equation can be approximated in the limit of r 3> 1 as 

P(r) « d'3(T-(^+2)+2er-(^+i))e-/-=(P). (17) 

In the scaling regime, i.e. er <C 1, we find that P{t) ^ 
T~" with a = (3 + 2. We performed the numerical sim- 
ulations with various values of /3 for both deterministic 
and stochastic cases. In Fig. [4] we find that the simula- 
tions confirm the analytic solutions for both deterministic 
and stochastic cases for the range of r < d. The range of 
T > d in the stochastic case is only qualitatively matched, 
which is again due to the approximations we have used. 

III. EMPIRICAL ANALYSIS 

The results of the models for different cases of time- 
varying priorities in Section |ll] can be applied to under- 
stand the empirical distributions of the updating times of 
papers that have appeared in arXiv.org (arXiv in short) 
and of the processing times of papers published in Phys- 
ical Review journals (PR in short). The essential ingre- 
dients of the empirical setup include the following: (i) 
An agent has a given task of interest, the execution of 
which gets postponed due to the other tasks and (ii) the 
priority of the task varies with time. 

A. Description of the datasets 

The first dataset consists of all the papers that ap- 
peared on arXiv between January 1995 and March 2010. 



For each paper we consider its updating time, defined as 
the time interval in days between submission and the fi- 
nal update. The papers that were not updated have been 
discarded in our analysis. This dataset contains 185151 
papers. We divide this dataset into three groups depend- 
ing on the number of authors in each paper, denoted by n. 
The groups of n = 1, n = 2, and n > 3 consist of 55382, 
56572, and 73197 papers, respectively. The papers with 
more authors have the smaller average updating time: 
220.8, 188.5, and 160.0 days for groups of n = 1, n = 2, 
and n > 3. 

The second dataset consists of papers that appeared 
in PR between January 2000 and December 2009. For 
each paper we consider its processing time, defined as 
the time interval in days between submission and publi- 
cation. We have retained only the research articles and 
have removed articles of other types such as errata and 
publisher's notes, which often have different processing 
times. The exact publication dates of papers printed be- 
fore 2000 were not available and hence were ignored. This 
dataset consists of 157484 papers. We divide it into two 
5-year periods, i.e. sets of papers appearing for 2000-2004 
and for 2005-2009, respectively. The average updating 
time decreased from 200.3 days for 2000-2004 to 163.9 
days for 2005-2009. We would like to note that many re- 
searchers submit their papers to arXiv at the same time 
while submitting them to journals and update the arXiv 
versions on acceptance of the article by the journal. Both 
the updating time and the processing time are denoted 
by T. 



B. Updating time distributions in arXiv.org 

We plot the updating time distributions of arXiv pa- 
pers for groups of n = 1, n = 2, and n > 3, as depicted in 
Fig.[5];a). In aU cases, the initial part of the distribution 
(t < 50) shows a power-law decaying behavior, which is 
followed by an exponentially decaying tail. The distribu- 
tion also shows a hump around r = 100 indicating its bi- 
modal nature. These distributions can be fitted by a sim- 
plified bimodal curve of the form P(r) oc {t~°' + c)e~'^^'^'' . 
We ignored very large updating times (r > 2000) due to 
the inconclusive fitting results. The parameter values 
are estimated as a « 0.76(7) and Tc ~ 234(29) for n = 1 
group, a sa 0.99(13) and Tc ~ 189(17) for n = 2 group, 
and a ~ 1.16(15) and ~ 166(7) for n > 3 group, re- 
spectively. The fitted curve for n — 2 group is plotted in 
Fig. [5]ja). It is found that the value of a is slightly larger 
for the papers with more authors. The cutoff updating 
times Tc are consistent with the corresponding average 
updating times within error bars. 

Let us consider the updating mechanism of the papers 
submitted to arXiv. Immediately after submission the 
authors are most alert and are likely to find some anoma- 
lies in their paper and subsequently update it. As time 
passes, they get more involved with other projects and 
the probability to find irregularities decreases, leading to 
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FIG. 5: (a) Probability distributions of updating times for 
papers by single authors {n — 1), two authors (n = 2), and 
more than two authors (n > 3) in arXiv.org. The solid line in- 
dicates the bimodal curve (r-°+c)e"^''^= fitting to the group 
of n = 2. The estimated values of parameters are provided 
in the text, (b) Waiting time distributions of the stochastic 
model with power-law decreasing priority in Subsection |II B| 
and its extended versions to n agents. Here r'"' is defined as 
min{Ti}i=i,... ,n with each drawn from the distribution of 
the single agent model. 



a decreasing priority of the updating task. In addition, 
there is a finite probability that the authors might receive 
a referee report, and hence will update the paper with a 
revised version even when the priority of the updating 
task is smaller than the other tasks. Thus, the decreas- 
ing priority and the stochastic selection of the updat- 
ing task can qualitatively describe the bimodal updating 
time distributions in arXiv. However, for the papers by 
more than one author the interaction among authors is 
also an important factor for determining the updating 
time. In this case, we consider the generalized model 
with n agents. One of the authors who first finds the 
anomalies in the paper is assumed to update the paper. 
Then, the resultant updating time is determined by the 
shortest updating time as r'"^ = min{Ti}i=i^... where 
each Ti is drawn from the distribution of the single agent 
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FIG. 6: (a) Probability distributions of processing times in 
Physical Review journals for different periods. We fit the in- 
creasing part and the tail part of the distributions with the 
curves (r" + c')e~^^'^'' and T~°'e~'^^^'' , respectively. The fit- 
ting curve for the tail of 2000-2004 is not shown for clear pre- 
sentation. The estimated values of parameters are provided 
in the text, (b) Waiting time distributions of the stochastic 
model with power-law increasing priority in Subsection |II C| 
and its extended versions to n agents. Here r'"' is defined as 
max{ri}i=i,... ,„ with each n drawn from the distribution of 
the single agent model. 



model in Subsection III bI The distribution of r*^"-' reads 

n 

P{t^"^) = c„P(ti = t(")) n Pin > t(")), (18) 

1=2 

where c„ is the normalization constant. We numerically 
obtain P(t("^) from the numerical results of P{Ti) of the 
single agent model for small values of n, as shown in 
Fig. [sj^b) . The decreasing average updating time accord- 
ing to n is observed as in the case of empirical results. 
The power-law exponent for the range of small updating 
times, denoted by a^m, tends to increase with n. This 
tendency is also consistent with the empirical observation 
apart from the exact values, while the discrepancy in the 
values of power-law exponent would require some addi- 
tional or different approach beyond the current version 
of our models. 
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C. Processing time distributions in Physical 
Review journals 

Next, we plot the processing time distributions of PR 
papers for periods 2000-2004 and 2005-2009, as depicted 
in Fig. [6]^ a). In both cases, when compared to the bi- 
modahty found in the arXiv case, the distributions are 
unimodal with peaks at t sa 150 and 120 for 2000-2004 
and 2005-2009, respectively. We fit the tail of distribu- 
tion, by using the functional form of P(t) c>c r~"e^'^/'^'=. 
We estimate a « 2.7(3) and « 234(27) for 2000-2004 
and a « 3.2(3) and w 337(38) for 2005-2009, respec- 
tively. While the values of a are the same within error 
bars, the value of Tc increases from 2000-2004 to 2005- 
2009. Then we fit the increasing part of distribution with 
the form P{t) oc (r" +c')e~'^/'^'' , where we have assumed 
the same value of Tc as for the tail. The estimated val- 
ues are a' « 3.2(2) for 2000-2004 and a' « 1.6(1) for 
2005-2009, respectively. The fitted curves are plotted in 
Fig.|6];a). 

The publication process of PR can be understood as 
a simple queuing process from the point of view of the 
referee. Among many factors affecting the processing of 
the paper, we consider only the time taken by the referee 
as the main factor. The other factors, such as the time 
taken by editors, are considered to be fixed and hence 
they are expected not to change the shape of the distri- 
bution but only to affect the location of the peak. The 
referee assigns some initial priority to the received pa- 
per, which may depend upon the quality of the paper 
and the referee as well. The priority of the reviewing 
task increases with time, as the journal editor hurries up 
the referee in case when he or she does not submit the 
review report within a stipulated period. The referees 
may also select the reviewing task randomly, indepen- 
dent of the priority. Hence the simple mechanism of the 
power-law increasing priority of the referee can describe 
the unimodal processing time distributions in the case of 
PR. As the number of referees per paper varies, we can 
consider the generalized model with n independent ref- 
erees. The processing time of a paper is determined by 
the longest reviewing time, i.e. r*^") = max{rj},-i ... „. 
Here each Ti is drawn from the distribution of the sin- 



following equation: 



+ 1 



gle agent model in Subsection |H C[ The distribution of 
P(r(")) reads 



P(t(")) - c„P(ti = t(")) J] ^(n < T 



(19) 



i=2 



where c„ is the normalization constant. When t'") ^ d, 
since P{Ti) ^ r" for each i, the leading term of P(r'-"-') 
becomes r'^"^""''^ with a'^ = na' -\-n — l. On the other 
hand, if t(") > rf, P(r, < r'")) is of the order of 1 due to 
the unimodality of the distribution. Thus, the tail of the 
distribution is described by the same exponent as in the 
single agent model, i.e. Omax — ot, which is independent 
of n. Then, by means of a = a' + 2 obtained in the 
unbounded case of the single agent model, we get the 



1 



(20) 



Since the information about the numbers of referees in 
PR is not accessible, we infer the effective number of 
referees by plugging the empirical values of power-law 
exponents a and a' into amax and Omaxi respectively. We 
obtain n « 2.5 for 2000-2004 and n « 1.2 for 2005-2009, 
respectively. The decreasing value of n is also consistent 
with the decreasing average processing time, in a sense 
that the longest reviewing time will be shorter in case 
of the fewer referees. The additional factors, such as the 
number of revisions needed and the time taken by the 
authors for revision, can be considered in an extended 
version of the model. 



IV. CONCLUSION 

Various forms of waiting time distributions have been 
observed in the human dynamics, including the bimodal 
and the unimodal distributions with power-law or expo- 
nentially decaying tails. To figure out the origin of vari- 
ous types of bursty correlations, we have studied a single 
agent queuing model with two tasks: one is the task of in- 
terest whose priority algebraically increases or decreases 
with time, and the other is the random task with the ran- 
dom priority. If we choose the model to be deterministic, 
only the highest priority task will be executed. In the 
stochastic version, however, one of the tasks in the queue 
is randomly selected for execution with small probabil- 
ity. In all the cases we obtain the approximate analytic 
solutions for the waiting time distributions of the task of 
interest, which are confirmed by the extensive numerical 
simulations. 

For the case of power-law decreasing priority with ex- 
ponent 7, we find that the power- law exponent of the 
waiting time distribution P{t) ~ t~" shows the non- 
monotonous behavior as 7 increases, i.e. a = 2 — 7 for 
7 <C 1 and a = 7 for 7 ^ 1. The stochastic case of the 
model leads to both the bimodality and the finite cut- 
off of distributions. For the case of power-law increasing 
priority with exponent /3, the resulting waiting time dis- 
tributions increase as power-law P{t) ~ up to some 
scale factor. The power-law exponent a of tails turns 
out to be 2 and (3 + 2 ior the bounded and unbounded 
case, respectively. The extended versions of our models 
to many agents are also compared to the empirical dis- 
tributions of updating times for the papers appeared in 
arXiv.org (arXiv) and of processing times for the papers 
published in the Physical Review journals (PR). The up- 
dating process in arXiv can be understood as a model 
where each author's priority of updating the paper de- 
creases with time as a power-law. The processing time of 
a paper in PR can be understood by means of the review- 
ing process, where each referee's priority of reviewing the 
paper increases with time. 
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Compared to the Barabasi's queuing model [S], in our 
models the tasks in the queue are assumed to be of dif- 
ferent types. The task of interest among them has the 
time-varying priority, while its execution is affected by 
the random tasks. Our simplified models for the human 
dynamics can be extended or generalized to the more 
complicated and realistic cases. As an example, a re- 
searcher working on two ongoing projects with different 
time-varying priorities can be studied by means of the 
model with two tasks of interest and one random task in 
his or her to-do list. 
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